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THE ROLE OF FROLOV’S CUBATURE FORMULA FOR FUNCTIONS 
WITH BOUNDED MIXED DERIVATIVE 

MARIO ULLRICH AND TINO ULLRICH 


Abstract. We prove upper bounds on the order of convergence of Frolov’s cubature formula 
for numerical integration in function spaces of dominating mixed smoothness on the unit 
cube with homogeneous boundary condition. More precisely, we study worst-case integration 
errors for Besov B[/g and Triebel-Lizorkin spaces F“g and our results treat the whole range of 
admissible parameters (s > 1/p). In particular, we obtain upper bounds for the difficult the 
case of small smoothness which is given for Triebel-Lizorkin spaces F[/g in case 1 < 6 < p < oo 
with 1/p < s < 1/6. The presented upper bounds on the worst-case error show a completely 
different behavior compared to “large” smoothness s > 1/8. In the latter case the presented 
upper bounds are optimal, i.e., they can not be improved by any other cubature formula. The 
optimality for “small” smoothness is open. 


1. Introduction 

The efficient integration of multivariate functions is a crucial task for the numerical treatment 
of many multi-parameter real-world problems. The computation of the integral can almost 
never be done analytically since often the available information of the signal or function / is 
highly incomplete or simply no closed-form solution exists. A cubature rule approximates the 
integral 1(f) = J^ 01 ^ d f(x)dx by computing a weighted sum of finitely many function values 
and the d-variate function / is assumed to belong to some (quasi-)normed function space 
Fd C U([0,l] d ) of continuous functions. The optimal worst-case error with respect to F^ is 
given by 

n 

h/)-X A ‘/U) - 

2—1 

In this paper we give an (almost) complete answer for the question on the correct asymptot¬ 
ical behavior of Int n (F^) for several classes of functions with dominating mixed smoothness. 
The by now classical research topic of numerically integrating d-variate functions with mixed 
smoothness properties goes back to the work of Korobov [25], Hlawka [22], and Bakhvalov [2] in 
the 1960s and was continued later by numerous authors including Frolov [12], Temlyakov [431 
05] 06] 09], Dubinin [7] [8j, Skriganov [39], Triebel [53], Hinrichs et al. pH] 17, EE], Hinrichs, 
Novak, M. Ullrich, Wozniakowski [18, T9;. Hinrichs, Novak, M. Ullrich [20], Dung, T. Ullrich 
[9], T. Ullrich [56], Dick and Pillichshammer [6], and Markhasin [28jES0nj to mention just a 
few. In contrast to the quadrature of univariate functions, where equidistant point grids lead 
to optimal formulas, the multivariate problem is much more involved. In fact, the choice of 
proper sets X n C [0, l] fi of integration nodes in the d-dimensional unit cube is the essence of 
“discrepancy theory” and connected with deep problems in number theory, already for d = 2. 
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(1.1) Int n (F^) := inf inf sup 

x 1 ,..., x n &[ 0 , l] d A 1 ,...,A ri eR ||/1|F rf < 1 
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We study Besov-Triebel-Lizorkin classes B® e and F® g with dominating mixed smoothness 
on the unit cube [0, l] d with homogeneous boundary condition and provide lower and upper 
bounds for Int n (A® () ) which are sharp in order. Here and in the sequel A stands for B or F. 
This represents the crucial step for getting the same order of convergence also for the larger 
periodic spaces A® 0 (T rf ) and non-periodic spaces A® g(R d ) as the recent paper m shows. As 
a motivation let us emphasize that Besov regularity is the correct framework when it comes 
to so-called kink functions, which often occur in mathematical finance; e.g., the pricing of a 
European call option, whose pay-off function possesses a kink at the strike price m ■ Indeed, 
the simple example f(t) = |2 1 — 1| belongs to B^ ^(T), but not to any of the Sobolev spaces 
W2(T), p > 1. Looking at the bounds (11.31) and El below we see that the main order 
of convergence only depends on the smoothness s. Hence, to explain the (optimal) order of 
convergence n ~ 2 for the numerical integration of the function /, the classes B® e represent an 
appropriate framework. The method of our proof allows us to tackle also the (technically more 
difficult) Triebel-Lizorkin spaces F® e , which are, in the special case 6 = 2, 1 < p < oo, the 
Sobolev spaces W®. Thereby, we obtain the first non-trivial upper bounds on Int n (W®) in the 
range p > 2, 1 /p < s < 1/2, if the dimension is larger than 2. 

We show that there exists a cubature formula, namely the classical Frolov construction, 
which provides the optimal order of convergence for “almost” the entire Besov-Triebel-Lizorkin 
scale of function spaces (including anisotropic mixed smoothness, see Remark 15.61 below). 
Note that in some cases we do not have matching lower bounds, cf. Theorem 11.21 however 
we conjecture the optimality also in those cases. Let us emphasize that, in contrast to well- 
known sparse grid approaches, the nodes and weights of this cubature formulas are constructed 
independently of the parameters of the respective spaces; i.e., independently of s, p and 6. In 
other words, we do not need to incorporate any further information on the function to be 
integrated. 

Frolov m introduced this method in 1976, see also [5lj for a tutorial paper on the method. 
For the definition of the method, let T £ M. clxd be a suitable matrix with unit determinant, see 
Section [2] and define the lattices X n = n^ 1 / rf T(Z d ), n > 1. Frolov’s cubature formula is then 
rather simply defined as 

(1-2) Qn(f) = - E /(*) 

xex n n[o,i) d 

and it is well known that |X n n [0, l) rf | x n. For the detailed definition of the cubature formula 
Q n and the spaces A® e we refer to Section [2] and Section 01 respectively. We define the 
worst-case error of the cubature rule Q n in a class F^ by 

e(Qn, F d ) := sup [ f(x)dx-Q n (f). 

|]/llF d <i 7[o,i]^ 

Furthermore, we define the spaces A® as the collection of all d- variate functions from A® „ 
with support in the unit cube, cf. (14.21) below. It seems that the cubature formula Q n of 
Frolov itself is suitable only for functions from A® 0 , see Remark 15.41 below. However, in 
a forthcoming paper, see EE we will present two general modifications of Q n that define 
cubature formulas with non-equal weights providing the same order of convergence also for 
periodic or non-periodic functions, respectively. Note that all upper bounds on Int n (A® 0 ) 
that will be proven in this article are constructive. In particular, we will prove all upper 
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bounds via Int n (F rj r) < e(Q n ,F ( i) for a specific cubature formula Q n . In fact, these cubature 
formulas will be Q n from (jl.2l) . for functions from A®g. 

The main results of this article are the following Theorems 11.11 & 11.21 


Theorem 1.1. Let 1 < p, 9 < oo (with p < oo in the F-case). 
If s > 1/p, then 


(1.3) 


Int n (Bpg) x n-^logn)^- 1 ^ 1 " 1 / 0 ) 


n > 2. 


If s > max{ 1/p, 1/6}, then 

(1.4) Int n (F^) x n-^logn)^ -1 ^ 1 - 1 /®) , n>2. 


The upper bound in relation (11.31) has been already proved by Dubinin [8], see also the list 
of references below. However, our method is different (and simpler) and gives a unified proof 
for Besov and Triebel-Lizorkin spaces. That is why we give also the relatively short proof for 
this case. 

The lower bounds for 9 = oo are due to Bakhvalov [3]. The technique in the recent paper 
Temlyakov [ 50] also gives the lower bound for 9 < oo. Our proof of the lower bounds in Section 
7 is based on atomic decomposition and similar to the approach in [9]. This approach also 
gives the correct lower bounds in case p < 1, see Theorem 17.31 

In addition, we particularly pay attention to the case of small smoothness which is present 
for the spaces F® g if p > 9 and 1/p < s <1/9. It has been already observed by Temlyakov [47, 
1351 4SJ that for the Fibonacci lattice rule on T J the asymptotical worst-case error order with 
respect to the Sobolev class W® differs essentially in the log-exponent from (11.31) in the critical 
range of parameters p > 2 and 1/p < s < 1/2. That is, Temlyakov [48] proved for this range 
that in dimension d = 2 


b n s (log bn) 1 s , if 1/p < s < 1/2, 

6“ s (log6 n ) 1_s (loglog6 n ) 1_s , if s = 1/2, 


where <f> n denotes the Fibonacci cubature rule with respect to the Fibonacci lattice with b n 
points, where b n denotes the ?ith Fibonacci number. A matching lower bound for arbitrary 
cubature formulas is conjectured to hold, but is still not proven. We were able to show a 
corresponding upper bound for F® e also in higher dimensions for the Frolov cubature rule. 


Theorem 1.2. Let 1 < p < oo and 1 < 9 < p < oo. 

(i) If 1/p < s < 1/9 then 

Int n (F^) < n-^logn)^" 1 )^- 5 ) , n>2. 

(ii) 11 s = 1/9 then 

Int n (Fp i0 ) < n _s (logn) (d_1)(1_s) (loglogn) 1_s , n>3. 

Concerning matching lower bounds there is so far only hope for e(Q n , F* e ). It is an interesting 
open problem to ask the same question for Int ra (Fp 6) ) in the critical range of parameters. In 
any case, we strongly conjecture that the given convergence rates are sharp. This would be 
another characteristic of a structural difference between the B and F-spaces, recently observed 
by Seeger and T. Ullrich [361 [3T|. in the mentioned parameter domain. 

Spaces of the above type have a long history in the former Soviet Union, see m m Eg S3 
and the references therein. The scale of spaces B® e contains two important special cases of 
spaces with mixed smoothness: the Holder-Zygmund spaces (p = 9 = oo) and the classical 
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Nikol’skij spaces (9 = oo). Note that Sobolev spaces with integrability 1 < p < oo, p / 2 
and s > 0 are not contained in the Besov scale. They represent special cases of Triebel-Lizorkin 
spaces if 9 = 2. The upper bounds in Theorem o were already obtained for 

P,(7 

• Wp(I 2 ), s > 1/p, see [37] Thm. IV.2.1 k IV.2.5] and the references therein, 

• B£ i0O (T 2 ), s > 1/p, see |47l Thm. IV.2.6], and the references therein, 

• B* 0 (T 2 ) for s > 1/p, see [BB] and [9] , 

• B/ fl ([0, l] 2 ) for 1/p < s < 1 + 1/p, see [56] and [9], 

• B* 0 ([O, l] d ) n {/: f(x ) = 0 if Xi = 1 for some i} for s < 1, see [28] , 

• F* 0 ([O, l] d ) (~l {/: f(x ) = 0 if Xi = 1 for some ?'} for max{l/p, 1/9} < s < 1, see [28] , 

• Wp([0, l] (i ) for s € N, see [T2], 55]. 

• Wp(T rf ) for 2 < p < oo and s > 1, see [37] Thm. IV.4.4], 

• B* !00 (T d ) if s > 1 and 1 < p < oo, see [37] Thm. IV.4.6], and 

. B;’ 0 (T rf )forl <p,9< oo and s > 1/p, see [7] [8]; 1/p < s < 2, see [17] . 

The last four results were obtained by using Frolov(-type) cubature rules. The other bounds 
are achieved for cubature formulas that use (digital) nets and the Fibonacci lattice rule, re¬ 
spectively, where the latter is restricted to d = 2. For the Besov spaces with d > 2 there are 
also some (not optimal) upper bounds by Triebel |53J for 1/p < s < 1 + 1/p using integration 
nodes from Smolyak grids or by Temlyakov [43] using quasi-Monte Carlo lattice rules of Ko¬ 
robov type. Additionally, some results for Triebel-Lizorkin spaces can be easily obtained by 
embedding into B* 0 , see Lemma 14.71 below. Concerning matching lower bounds for Sobolev 
spaces W* we refer to [45] [48], and for Besov spaces B^ g in case s > 1/p and 1 < p, 9 < oo 
to the recent paper [9]. This will be complemented in Section [7] for the quasi-Banach case 
min{p, 9} < 1. It shows that there will be no additional gain in the convergence rate when 
passing to the situation p < 1, an observation particularly relevant for the above mentioned 
kink functions. In fact, we will present asymptotically optimal results for 

1) Int n (A* 0 ) in the quasi-Banach cases min{p, 9} < 1, see Section [777T1 and 

2) Int n ( A .^//) in the limiting case (s = 1/p), see Section [6721 

Note that the results in Section 16.21 require additional assumptions on 9 to assure continuity 
of the functions. For the precise statement consider the mentioned section. 

In addition, it is worth mentioning, that only recently Krieg and Novak [26] designed a Monte 
Carlo algorithm that is based on the Frolov cubature rule and proved that this algorithm gains 
the order 1/2 in the main rate compared to the deterministic version studied in this paper. 

The paper is organized as follows. In Section [2] we introduce the Frolov cubature rule 
in detail. In addition, we give a direct construction of the Frolov lattice matrix without 
matrix inversion. Afterwards, in Section [3] we collect several tools from harmonic analysis. 
In particular, Poisson’s summation formula and Calderon’s reproducing formula in connection 
with local mean characterizations of function spaces turn out to be crucial for the error analysis 
of the cubature formula. The function spaces of interest will be introduced in Section [4] In 
Sections [7] and [B] we prove the upper bounds for Frolov’s cubature rule in the spaces A* d for 
the Banach and the quasi-Banach situation, respectively. The matching (except for the case 
of small smoothness) lower bounds will be proven in Section 0 

Notation. As usual N denotes the natural numbers, No = N U {0}, Z denotes the integers, 
R the real numbers, and C the complex numbers. By T we denote the torus represented by 
the interval [0,1], where the end points are identified. The letter d is always reserved for the 
underlying dimension in R^Z^,!^ etc. We denote with ( x,y) or xy the usual Euclidean inner 
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product in W l and C d . For a £ M we denote a + := max{a,0} and a = |_aj + {a} with [oj £ Z 
and 0 < {a} < 1. For 0 < p < oo and x £ we denote \x\ p = (Xa=i |»Ci| p ) 1/,p with the usual 
modification in the case p = oo. By (xi ,..., Xd ) > 0 we mean that each coordinate is positive 
and [xj (resp. {x}) are meant component-wise. By T we denote the torus represented by the 
interval [0,1]. If X and Y are two (quasi-)normed spaces, the (quasi-)norm of an element x 
in X will be denoted by ||x||jv- The symbol X Y indicates that the identity operator is 
continuous. For two sequences of real numbers a n and b n we will write a n < b n if there exists 
a constant c > 0 such that a n < cb n for all n. We will write a n x b n if a n < b n and b n <a n . 


2. Frolov’s construction 


In this section we introduce the matrix T that is used in the construction of the lattice X n for 
the cubature rule Q n from (11.21) . We follow to a large extend the original work of Frolov 
see also mM - Let 

d d 

( 2 . 1 ) P d (t) := l[{t - 2j + 1) - 1 = ]J(t - £j), t e R, 


3 = 1 


3 = 1 


for suitable i.e. the are the roots of P d . It is easy to see that the are indeed 

real numbers and all different. Now define the matrix 

(2-2) f = 

and let T n = (ndet(T)) T such that det(T n ) = 1/n. The lattice for the cubature formula 
Qn is given by X n = T n (1 d ). 

Later we will see that the error of the cubature formula Q n applied to some function / is 
given by a sum over function evaluation of the Fourier transform Tf at the points B n ( Z d ), 
where B n is given by 

(2.3) B n = (T- 1 ) T . 

The lattice B n { Z rf ) is usually called the dual lattice of X n . The following result contains the 
most important property of B n . 

Lemma 2.1. Let B n be as above. Then, for each z £ B n (X d ) \ {0}, we have 


(2.4) 


d 

n 

3 = 1 


z i\ Z n- 


Proof. First note that it is known that inJ=i(Tm);| > 1 for each m £ \ {0} with T from 

m , see m Lemma IV.4.3] or [5T . Lemma 8]. Since T is a Vandermonde matrix we can 
write T~ l as the product where H is a Hankel matrix with integer entries and 

det(H) = —1, and D is a diagonal matrix with det(H) = det(T) 2 , see [27, Section 4]. Hence, 
d d 


inf 

m£Z d \{ 0 } 


n |(S n m) 7 -| = n det(T) inf TT dD^TH^m) 

J /n't rr 


i=l 


m£7j d \{ 0} 


n 


inf 


> 


det(T) m&L d \{ 0} 
n 


IJ\(TH t v^\ > 


n 


inf 


det(T) m£Z d \{ o} f ^ 


n im 


det(T) 
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where we have used that H T (7L d ) C Z rf . □ 

Note that lattices which satisfy the conclusion of Lemma 12.11 satisfy |X n D [0,1)°| = n + 
0(ln d_1 (n)), see [39] , 

Remark 2.2. The original construction of the lattice X n of Frolov [12], see also m, uses the 
generator n~ l ^ d B\ instead of T n . As it turns out we can use (|2.2I) directly and do not have to 
invert the matrix T for the construction of the nodes. However, this “trick” only works if the 
generator T is (a multiple of) a Vandermonde matrix. 


The crucial property of the matrix B n , which we will need in the sequel, is reflected best by 
the numbers 


(2.5) 


Z n (m) 


(B n (z d ) \ {o}) n i m 


with 


(2.6) I m := {x £R d : C\ < \xj\ < C 2 2 mj for j = 1,... , d} , m £ Nq, 

with constants 0 < C± < C- 2 < oo independent of m. The numbers Z n {m ) denotes the number 
of points from B n (Z d ) excluding 0 in (the union of at most 2 d ) axis-parallel boxes with volume 
of approximately 2l m l 1 . The following estimates for Z n {m ) are a direct consequence of (12.41) 
and the lattice structure of B n {Z d ). For a proof we refer to [37[ IV.4] or [53] . 


Lemma 2.3. There exists an absolute constant c < oo such that with 


(2.7) r n := log 2 (n) — c 

the numbers Z n (m) from (12.51) . n £ N, mn £ Nq, satisfy 

(i) Z n (m) = 0, if \m\i < r n , and 

(ii) Z n (m ) < 2l m l 1 /n, otherwise. 

Remark 2.4. The choice of the polynomial (12.11) is quite flexible. One could replace by 
every irreducible (over Q) polynomial with integer coefficients that has d different real roots. 
Another example, for d = 2 e , is the Chebychev-type polynomial P 2 ?.(t) = 2 cos(2^ arccos(t/2)), 
see mi p, 242]. The roots of these polynomials are given by = 2cos(7r(2i — l)/2 n+1 ), 
i = 1,..., d. 


Remark 2.5. It is worth noting that the results of this paper do not rely on the specific 
construction. In fact, every lattice for which the dual lattice satisfies the bound of Lemma [2. II 
would lead to the same results. See [39] for a detailed study of such lattices. 


3. Tools from Fourier analysis 

3.1. Preliminaries. Let L p = L p (R d ), 0 < p < oo, be the space of all measurable functions 
/ : —> C such that 

\\f\\p '■= ( / \f(x)\ p dx) /P < oo 
v JR d J 

with the usual modification if p = oo. In addition, we denote by C(M d ) the space of all bounded 
and continuous comp lex-valued functions on M. d . 

We will also need L p -spaces on compact domains 9 C instead of W l . We write ||/||l (n) 
for the corresponding (restricted) L p -norm. For / £ L i(M rf ) we define the Fourier transform 

FI® = [ f(y)e- 2 ^ y) dy, ^R d , 

JR d 
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and the corresponding inverse Fourier transform .F _1 /(£) = Ff (—£). Additionally, we define 
the spaces of continuous functions C(R d ), infinitely differentiable functions C°°(M. d ) and in¬ 
finitely differentiable functions with compact support Cg^R^) as well as the Schwartz space 
S = <S(R rf ) of all rapidly decaying infinitely differentiable functions on R d , i.e., 

<S :={</? £ C°°(R d ): < oo for all k,l £ N} , 


where 


IMIm 


a + M)‘ E i^ (m) (-)i 

m =0 


oo 


The space S'(M. d ), the topological dual of 5(R d ), is also referred to as the set of tempered 
distributions on R d . Indeed, a linear mapping / : <S(R d ) —> C belongs to S'(W l ) if and only if 
there exist numbers k,i £ N and a constant c = Cf such that 


(3.1) |/(v>)| < Cf\\(p\\ k ,e 

for all p £ 5(R rf ). The space <S'(R rf ) is equipped with the weak*-topology. The convolution 
p * ip of two square-integrable functions ip, ip is defined via the integral 

x)= p(x - y)ip(y) dy . 

JR d 

If p, ip £ 5(R rf ) then p * ip still belongs to S(W l ). In fact, the convolution operator can be 
extended to <S(R d ) x L\ , in which case we have p * ip £ <S(R rf ), and to 5(R rf ) x ^(R^) via 
(p * f)(x) = f(p(x — •)). It makes sense point-wise and is a C^-function in R d . As usual, the 
Fourier transform can be extended to S'(R d ) by ( Ff)(p ) := /(Tp), where / £ ^(R^) and 
p £ <S(R rf ). The mapping F : 5'(R d ) —> ^'(R^) is a bijection. 


3.2. Periodization and Poisson’s summation formula. The analysis of the error of cu- 
bature formulas that use nodes from a lattice is naturally related to an application of Poisson’s 
summation formula and variations thereof, see (15.211 . A more detailed treatment and a proof 
of the following theorem can be found; e.g., in [4Qi Thm. VII.2.4 & Cor. VII.2.6]. 

Theorem 3.1. Let / £ Li(R). Then its periodization /(^ + •) converges in the norm of 

Li([0, l] rf ). The resulting (1-periodic) function in Li([(), l] rf ) has the formal Fourier expansion 

Ff(k ) e 2mkx . 

keZ, d 

Moreover, if / £ S( R d ), then both sums converge absolutely, and hence point-wise. 

Note that absolute convergence in the above sums actually holds under weaker assumptions 
on /. However, this version is enough for our purposes. 

In what follows, cf. (15.21) . we will need a point-wise version of Poisson’s summation formula, 
like Theorem 13.11 in cases where we cannot guarantee that the sums converge absolutely. In 
such cases we have to specify in which sense (ordering) we understand convergence. This is 
provided by the corollary below. 

Corollary 3.2. Let / £ Li(R d ) be continuous with compact support, T £ W dxd be an in¬ 
vertible matrix, and B = (T _1 ) t . Furthermore, let po £ Cg°(R) with <^q( 0) = 1 and de¬ 
fine pj(t) := p 0 (2~ j t) - p 0 (2~ j+1 t), j £ N, t £ R, as well as the (tensorized) functions 
Pm(x) := p mi (xi) • • • • • Pm d (x d ), m £ N(j, x £ R d . Then 

det(T) V f(T£) = lim V V p m (Bk) Ff (Bk). 

N^roo 

££7, d m: Imloo^W k£Z d 
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In particular, the limit on the right hand side exists. 

Proof. The proof is based on [40l Theorem VII.2.11], We put <I > _b(-) := po(B-) and note, that 
satisfies the assumptions [301 (2-10)], in particular $_b( 0) = 1. Moreover, / is continu¬ 
ous and has compact support, which implies that the periodization det(T) Yle£Z d f{T(£ + x)) 
belongs to (^(T^). Applying |4H1 Thm. VII.2.11] together with Theorem 13.11 above we obtain 

det(T) V f(T{£ + x)) = lim V $ s (2 ~ N k) Ff (Bk) e 2nikx 

z — 4 N^oo z — 4 

l&A k£L d 

in C(T fZ ). Setting x = 0 gives 

(3.2) det(T) y f(T£) = lim V <h B (2 ~ N k)Tf(Bk). 

z —' N—too z — 4 

By construction we note that 

^b{^~ N k) = ipo(2~ N (Bk)i) • ... • ipo(2~ N (Bk)d) = E Pm(Bk). 

m: |m|oo<AI 

Plugging this into (13.2D and interchanging the order of summation yields the desired result. □ 


To prove our main results we will also have to bound the norm of certain series of functions. 
In fact, we treat two different kinds of functions: The first bound in Lemma 13.31 requires that 
the functions itself have compact support, while the second one in Lemma 13.51 is for functions 
with compactly supported Fourier transform. 


Lemma 3.3. Let B G M dxd be an invertible matrix and G be a bounded set. Further¬ 
more, let {/ m } meN d C <S(R rf ) be functions with supp(/ m ) C 17 for all m G Nq and define 


M B , n ■= {l&Z d :(£+ [0, l) d ) D B T (Q) + 0} 
Then, for 1 < 8,p < oo, we have 


E | 

mgHj 

In particular 


2ni £• 


meNg 


6K 1/0 


y/(w)e : 


2tt i£. 


< 


< 


Mr o \ 1 ~ 1 /p 


det(-B) 


E 


meNg 


i/e 


Mb a \ 1 ~ 1 /p 


det (B) 


v 


for / G 


with supp(/) C Q. 


Proof. Let p' be given by 1/p' = 1 — 1/p. Let h m (x) = f m (Tx) with T = (I? _1 ) T and note 
that supp (h m ) C B T ( 12). Clearly, J~ f m {Bt) = det (T)/Fh m (l) and, hence, by Theorem 13.II and 
Holder’s inequality we obtain 


E |E^W eW ' 


0 \ i/e 


jtiSNq t&A 


L P ([0,l] d ) 


det ( T ) E Em*+-: 


< det(T) E ( E \hm(£ + 


iez d me Ng 


6\ i/e 


< det(T) 


M 


W 

B,n 


E ( E I hm(£+-) 


meNg l&A 


L P m\ d ) 

9\ p/9\ 1/p 


6 \ 1/0 


ip([0,l] d ) 


i&1 d meNg 

Performing the integration and interchanging sum and integral yields 


MIOdP) 
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E I E Ffmm' 


2ni£- 


mGNg feZ d 


6K 1/6 


ip([0,l] d ) 




E (E 

e&Z d mGNg 


(E / 


E 

h“m (*C 



+ [0,1]1 

\ * ^ 

trENq 


<° r (L 


hm (^) 

e^p/e 



me Ng 



1M b , 

\ i-i Ip 

n J 

(E \u e ' 

.i/e 


d\p/e ^i/p 

dx 


9\p/0 \ i/p 

ax 


dx 


The second statement follows if we set fo = f and f m = 0, m/ 0. 


□ 


Remark 3.4. Note that Mbq is the number of unit cubes in the standard tessellation of R d 
that are necessary to cover the set H T (fi), while det(R) equals the volume of JJ T ([0, l] rf ). This 
shows that with a matrix of the form B n = rd^B, n > 1, cf. (12.31) . we obtain 


lim 

n—>oo 


M B n ,n 

det (B n ) 


for every Jordan measurable set fl. 


vold(fi) 


Lemma 3.5. Let B G M. dxd be an invertible matrix and fl G R rf be a bounded set. Further¬ 
more, let g G <S(R rf ) with supp^g) C fh Then, for 1 < p < oo, we have 


E F9(Bk)e 


2nik- 


k^Jj d 


L P (m d ) 


< 


b( z d ) n n 


i-i /p 


Ibll 


Proof. Clearly, the proof of Lemma 13.31 for p = 1 works also for functions f m without compact 
support and, hence, already proves the statement for p = 1. (Set fo = 9 and f m = 0 for 
m / 0.) For p = oo we easily obtain the upper bound \B(Z d ) nfi| • ||.7 : ’g'|| 00 < \B(Z d ) nfi| • ||g||i. 
Hence, the lemma follows by Holder’s inequality. □ 


3.3. A discrete version of Calderon’s reproducing formula. Our analysis heavily relies 
on a discrete version of Calderon’s reproducing formula [5], eq. (3.1)]. A “continuous” and 
homogeneous version of the Lemma below has been proved in [23]. This principle has been 
used by several authors 0, mi, is], m to prove equivalent (local mean) characterizations 
for Besov-TViebel-Lizorkin spaces, see Section [4] 

Lemma 3.6. Let To, 1 ! 7 ! G <S(R) be functions with 


(3.3) 

|^o(0l>0 for |£| < e 

and 


(3.4) 

l^!(OI >0 for | < |e| < 2e 


for some e > 0. Then there exist Aq, Ai G <S(R) such that 
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(i) supp-FAg C {t G M : \t,\ < e} 

(ii) suppJ-'Ai C {t G M : e/2 < |t| < 2e} and 
(Hi) for all £ G M, 

OO 

(3.5) = 1, 

3=0 

where T.,(x) = 2- J ~ 1 T] (2 J ~ l x) and A j(x) = 2- ?_1 Ai(2 :,_1 x) for j G N. 


Proof. Following [58| Thm. 1.20] we use the special dyadic decomposition of unity with ip(t) = 1 
if \t\ < 4/3 and ip(t) = 0 if \t\ > 3/2. Put To := and 3>i := 2To(2-) — To, i.e. = 

Tq(-/ 2) — T 0 . With Tj := 2 J ' -1 Ti(2- 7 ' _1 -) for j > 1 we define Aq, Ai through 


TAj(t) := 


F<j>j(2t/e) 


t G M. 


□ 


We dehne the d-fold tensorized functions 

d d 

A m (x) ■■= JjA mi (xj) and T m (x) := JJ T mi (xj) , x£R d . 

i =1 i=l 

where A j, 4 !j, j G N, are defined in Lemma 13.61 We obtain from (13.511 the identity 
(3.6) £^(0-^(0 = ! , 

meNg 

By the construction of the tensorized functions (and Lemma [3.611 we know that the support of 
J-A m is of the form (12.61) and we will write in the sequel 

I m := supp J r A m , meNj}. 


4. Function spaces with dominating mixed smoothness 


In this section we introduce the function spaces under consideration, namely, the Besov and 
Triebel-Lizorkin spaces of dominating mixed smoothness. Note that the Sobolev spaces of 
mixed smoothness appear as a special case of the Triebel-Lizorkin spaces. There are several 
equivalent characterizations of these spaces, see [58] . For our purposes, the most suitable is 
the characterization by local means (see [ 581 Theorem 1.23] or [55i Definition 2.5]). 

We start with the definition of the spaces on R d . 

Let To, Ti G Co°(M) be such that 

(*) \^M0\ > 0 for ICI < 

(ii) |J r Ti(^)| >0 for | < |£| < 2e and 
(in) i(0) = 0 for all 0 < ct < L 

for some e > 0. A suitable L will be chosen in Definitions 14.21 & m As usual, we dehne 

Tj(x) = 2 j ~ 1 ^ 1 (2 j ~ 1 x), j G N, 


and the (d-fold) tensorization 
(4.1) 


d 

T m (x) = JJT mi (Xi), 
i= 1 


where rn = (mi,..., rn f i) G Nq and x = (x \,..., Xd) G R d . 
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Remark 4.1. There exist compactly supported functions To,'hi satisfying (i)-(iii) above. 
Consider To to be the up-function, see Rvachev [33] . This function satisfies To € Co°(M) with 
snpp(To) = [-1,1] and .FTo(£) = n^Lo sinc(2 k £), £ £ R, where sine denotes the normalized 
sinus cardinalis sinc(£) = sin(7r£)/(7r£). If we define Ti € Cq°(M) to be 

d L 

*i(®) := d^( 2 T°( 2 ‘) -T°(0)(*), 

it follows that .FTi(£) = (2m£) L (j r To(^/2)—J r To(^)). It is easily checked that these functions 
satisfy the conditions (i),(ii),(iii) above. In particular, (i) and (ii) is satisfied with e = 1. 
Moreover, we have for all mGNjJ that the tensorized functions T m satisfy supp(T m ) C [—1, l] d . 
We will work with this choice in the sequel. 

Let us continue with the definition of the function spaces A® g = A s p g (U. d ) with A € {B.F} 
defined on the entire M. d . 

Definition 4.2 (Besov space). Let 0 < p, 9 < oo, s E M, and {T m } mgN d be as above with 
L + 1 > s. The Besov space of dominating mixed smoothness B® e = B® g (R d ) is the set of all 
/ € <S'(M d ) such that 

II/IIbj,. := ( £ 2*Mi«||* m »/||») 1/ ' < oo 

meNg 

with the usual modification for 6 = oo. 


Definition 4.3 (Triebel-Lizorkin space). Let 0 < p < oo, 0 < 9 < oo, s € M, and {T TO } mgN d 
be as above with L + 1 > s. The Triebel-Lizorkin space of dominating mixed smoothness 


F® 0 = F® e (M d ) is the set of all / 6 <S'(M d ) such that 


■ p,« 


( E 2 s|m|lfil |T m * /(•)| 0 ) 1/6 
me 


< OO 


with the usual modification for 9 = oo. 


Remark 4.4. In the special case 9 = 2 and 1 < p < oo we put W® := F® 2 which denotes the 
Sobolev spaces of dominating mixed smoothness. It is well-known (cf. |38l Chapt. 2] or [58] ), 
that in case s £ No the spaces W® can be equivalently normed by 

ll/llwj - ( £ l|o“/|g) 1/p . 

aeNg 

|a|oo<s 

Remark 4.5. Different choices of To,Ti lead to equivalent (quasi-)norms. In fact, it is not 
even necessary that To and Ti have compact support. However, for the proof of our results 
this specific choice is crucial. 

Remark 4.6. Note, that the spaces B® e (R d ) and F® d (R d ) are usually defined via a dyadic 
decomposition of unity on the Fourier side (which represents a special case of the above frame¬ 
work) like the one we used in the proof of Lemma 13.61 see [58] . Here we used a different 
characterization which was proven to be equivalent, see |58j and [55] . Let us comment on the 
recent and non-trivial history of those characterizations. 

• In 1992 Triebel proved those characterizations for the isotropic spaces F^ q (R d ), see [52l 
2.4.2, 2.5.1]. He obtained characterizations for p, 9 > 1. 
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• Later, in 1999, Rychkov [34] extended it to the quasi-Banach case. However, there is a 
gap in his proof (observed 2007 by Hansen), see [14| and Remark 4.4 in [55] . 

• In 2005 Vybiral modified Rychkov’s method for the dominating mixed case. Vybfral’s 
proof also contains Rychkov’s gap. However, Vybiral was the first who did the local 
means for dominating mixed spaces. 

• The proof in [55] differs from Vybfral’s proof and fixes the mentioned gap, see also 
Hansen HU. The proof is based on a maximal function technique, see [55., pp. 20], due 
to Stromberg, Torchinsky m Chapt. V] which has been already proposed by Rychkov 
[351 Thm. 3.2], 

• The proof in [55] is a bit more general, namely for function spaces on semi-infinite 
rectangular domains. 

The next lemma collects some frequently used embedding properties of the spaces. 

Lemma 4.7. Let 0 < p < u < oo (p, u < oo in the F-case), s, t £ R, 0 < 9, p < oo. 

(i) For equal p we have the chain of embeddings 

B p,min{p,e»} ^ F p,e ^ B p,max{p,6>} 

(ii) In addition, whenever s — 1/p = t — 1/u the following “diagonal embeddings” hold true 

B p ,0 ^ B u,6» j F p,e ^ 

(in) as well as the embeddings 

F ;,e^Kp > 

Proof. For a proof we refer to [38, Chapt. 2], as well as [43] and [15] . □ 

Remark 4.8. The embeddings in Lemma l4.7l iii) are highly non-trivial. Let us comment on 

the history of these important embeddings. Although they are referred to as Jawerth-Franke 

embeddings several different people contributed to this result, especially what concerns Besov- 

Lizorkin-Triebel spaces with dominating mixed smoothtness. 

• An analog of the first embedding in Lemma [4.7H ii) for isotropic Besov-Triebel-lizorkin 
spaces Fp q (R d ) and B s pq (\R d ) has been obtained by Jawerth [24], the second one by 
Franke [10] . The proof was based on real interpolation. 

• A new proof of both relations based on atomic decompositions has been given recently 
by Vybiral [59] without using real interpolation techniques. 

• Hansen and Vybiral jT5] extended this technique to prove the respective embeddings 
for spaces with dominating mixed smoothness B® 0 and F® g . Note, that the technique 
used by Jawerth and Franke (using real interpolation) does not apply for spaces of 
mixed smoothness. 

• A special case of the embeddings in Lemma l4.7l iiif with the appropriate Sobolev type 
spaces W®(T d ) (where 9 = 2) instead of F® g is due to Temlyakov |42, ii], 

• The second embedding in Lemma I4.7f iii) in the special case when F® g is replaced by 
L u in the univariate univariate case goes back to Ul’yanov m and Tinian m- 

In the sequel we will always assume that s > 1/p This assures that the functions in B® g , 

F® e and W®, respectively, are continuous, see [38* Chapt. 2], With the same reasoning we 

obtain that B l J g -A C{R d ) for 6 < 1, and F X j v g -A C(R d ) for p < 1. 
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The class of functions we are interested in throughout this article are the subclasses of 
functions from A^ g (M d ) which are supported in the unit cube [0, l] rf , i.e. we consider the 
classes 

(4.2) A* Pi0 := {/ 6 A; i0 (M d ): supp(/) C [0, l] d } 

for A € {B, F}. 

5. Integration of functions from A p0 

Here, we study the cubature formula Q n from (11.21) . see Section [2], for the spaces A s e defined 
in (I4.2p . Our main results read as follows. 

Theorem 5.1. For each 1 < p, 9 < oo and s > 1/p, we have 

e(Qn, ty,o) X n _s (log 7j,)( d-1 )( 1-1 / 0 ). 


Theorem 5.2. For each 1 < p < oo, 1 < 6 < oo and s > 1/p the following bounds hold true. 

(i) If s > 1/d then 

e(Qn, F s p ,e) - n-'aogn)^- 1 )^- 1 ^. 

(ii) If s < 1/9 then 

e(Qn, Flo) $ n-‘Qogn^-W-'K 

(Hi) If s = 1/9 then 

e {Qn,~F s P) e) ^ n _s (logn) (d_1)(1 “ s) (loglogn) 1_s . 


Setting 9 = 2 immediately implies the following important special case. 
Corollary 5.3. For each 1 < p < oo and s > 1/p the following bounds hold true. 

(i) If s > 1/2 then 

e(Qn,Wp) x n~ s (logn) (d ~ 1)/2 . 

(ii) If s < 1/2 then 

e(Qn, W-) < n-{logn)^ 1 * 1 -'). 

(in) If s = 1/2 then 

e{Qn, Wp) < n~ s (log n) ( ^ 1)/2 Vlog log n . 


Remark 5.4. The subclasses A* e of functions with support inside the unit cube [0, l] d rep¬ 
resent the “natural domain” for the Frolov cubature rule. This comes from the fact that the 
proof heavily relies on the use of Poisson’s summation formula (Corollary 13.2p . which requires 
that we can replace the summation in (11.21) over X n n [0, l) d by a summation over X n without 
changing the value of the sum, see m ■ 
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Remark 5.5. Let us give some historical comments on the special case in Corollary 15.31 
This result was partially known before and has a long history, starting with the work of 
Bakhvalov |2j. He proved a version of Corollary I5.3f ii for the Fibonacci cubature rule d> n 
instead of Q n for the periodic spaces W*(T 2 ) in case s£N and d = 2. This was extended by 
Temlyakov to 1 < p < oo. In fact, Temlyakov [461 [48] proved a version of Corollary I5.3f i)— (iii) 
with Q n replaced by for periodic spaces W*(T 2 ) in case d = 2 and s > 1/p. Corollary 15.31 7) 
in case d > 2, p = 2, and s € N is due to Frolov [12] which trivially implies the correct order for 
p > 2 due to embedding. For d>2, l<p<oo and s £ N, we refer to Skriganov [39]. All the 
remaining cases for d > 2 are new. In addition, the results for small smoothness, see (ii) and 
(iii), have not been known before. Also the recent paper Hinrichs, Markhasin, Oettershagen, 
T. Ullrich [T7] provides the correct order of Int„,(Wp(T rf )) in case 1 < p < 2 and 1/p < s <2 
with a quasi-Monte Carlo method based on order-2 digital nets. For more detailed historical 
comments see the monograph/survey Temlyakov [47149 ]. 

5.1. The proofs. We begin this section with the derivation of an error formula that is the 
starting point for the proofs in the specific cases. This explicit formula for 1 1(f) — Q n (f) I 
follows immediately from the point-wise version of Poisson’s summation formula for general 
lattices, see Corollary 13.21 

In this section we consider / € A* e , i.e. functions with support in the unit cube. Hence, we 
can rewrite our cubature rule Q n from (1 1.2 [i as 

(5-1) Q n {f ) = - 22 /( x ) = det ( T «) E 

aigx n n[o,i) d zex„ 

In the sequel we will use the specific kernel from Remark 14. II and its tensorized version T m , 
m € Nq . The corresponding functions A m , m € Ng, are given by Calderon’s reproducing 
formula (13.51) . (13.61) . The construction of these functions, cf. the proof of Lemma 13.61 assures 
that the functions p m := J[A m * ^> m ] satisfy the assumptions of Corollary 13.21 and hence, 

Qn{f) = 22 22 ^m*^m](B n k)Tf(B n k). 

mgNg A;SZ d 

Note that the inner sum is finite and that, actually, the outer sum is defined as a certain limit, 
see Corollary 13.21 However, we will see that this sum converges absolutely in all cases under 
consideration. 

Note that (e 2mk ', e 2mi '} = 1, if k = £, and 0 otherwise, where (•, ■) is usual inner product in 
L 2 ([0, l] d ). Using this together with 1(f) = J 7 f( 0) we obtain 

W)~Qn(f)\ = | XI Y. :F ^m{ B nk)^ m (Bnk)Ff(B n k) 

MO 

(5.2) = | E EE FK m (B n k)FYH m * f](B n t) (e 2 ^',e w -)| 

mgNg MO MO 

= I E (E ( s **o e2nik ' w) ) • 

mgNg MO MO 

Now we have to proceed differently depending on the space under consideration. In fact, we 
will perform in either case Holder’s inequality twice, but in a different order. 
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The result for B* fj . We now prove Theorem 15.11 
Using (15.211 we obtain by Holder’s inequality 


meNg IjtO 


2ni£- 


< E IIE- 77 ^^) 1 


27 rik- 


mSNg MO 


V([°d] d ) 


Y,n^m*f](B n l)e^ ie - 




Lp([0,l] d ) 


with 1/p + 1/p' = 1. 

Using Lemma 13.51 for the first and Lemma 13.31 for the second factor, we obtain 

W)-QnU)l< £ Znim) 1 - 1 '"'y 1 '” IIA^IItII*",»/II, 

me Ng ^ n 

By construction, the third factor is bounded by a constant. The second factor converges, as 
n —y oo, to vold(supp(\b m * /)) 1_1 /p < vol^Q—1,2] d ) 1 “ 1 / p , see Remarks 13.41 fc ITTl Hence, we 
obtain with Lemma 12.31 that 

(5.3) __ 

W)~Qn(f)\ < E (2 |m|l /n) 1/p ||tf m */|| P = n- 1 /? E 2 Hi(i/p-) 2 -Hi||^ m + /|| p 

mSNg: mGNg: 

|m|i>r n |m|i>r n 

with r n = log 2 (n) — c from (|2.7D . Applying Holder’s inequality one more time, with 1/6 + 1/6' = 
1, we finally obtain 


e(Q n , B s p>0 ) < E 2 e '^ 1 ^y /9 ' < n "V* 

m: |m|i>r n 

< n-^logn)^ 1 ^ 1 - 1 / 0 ), 


E (i+i) d - i 2 e ' e - { * 1 /p~ s )y /e 

i>r n 


since s > 1/p. This proves Theorem 15.11 


The result for () . We prove Theorem 15.21 If s > max) I jp, 1/9} Theorem 15.21 fil directly 
follows from Theorem l5.ll In fact, in case p <6 we have g B* e for all s. If p > 6 we use 
the embedding F* e s due to the compact support and the definition via local means. 
To apply the results in Theorem 15.11 we need s > 1/9. It remains to deal with the situation 
p > 6 and 1/p < s <1/6. 

Again, let 1/p + 1/p' = 1/6 + 1/6' = 1. We obtain from (15.21) that 


\I(f)~Qn(f)\ = | E f](B n £)e 

£^0 


2n ii- 


mgNg M 0 


£ 


m£Ng 


'[0,l] d 


E -FA m(B n k ) e 2 ^ ( E * /] (BJ) e“) d® 

M o e+o 
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Interchanging summation and integration and applying Holder’s inequality twice yields 

e\ 1/0' 


|/(/)-Qn(/)|< / ( V 2-^H 1 ! YFK m {B n k)e 2mkx6 ') 

■W MO ^ 


E 2 ®‘ Hl | E e“|') 1,9 dx 

neNj( MO 

e\ i/0' 


< 


E 2~ e ' s ^\Y^ m (B n k)e 

meNg MO 


27 rikx 


V([°>i] d ) 


E 2 0s Hi| Y^m*f](Bn^e 


2ni£x 


mgNg 


0 \ 1/0 


M([o,i] d ) 


At first, we bound the second factor. For this recall that the supports of *1 > m * / are subsets 
of [— 1,2] rf and set M n = M Bn ^_ _ 12 ]d in Lemma 13.31 Hence, we obtain from Lemma 13.31 (with 
fm = 2 s Hi\I/ m * /) that 

E 20SHl I E * /K 5 ^) ^ ^ 1/9 


meNg 


M([o,i] d ) 


< (det(T n )Af n ) 1 1/p ||( E 2 ' SHl 


* m*f 


e \ l / d 


meNg 


< 

r^j 


P ,e 


J2nikx 


It remains to bound the first factor. For this let 

A m ,n(x) := E^mC-Bn^je 2 

M 0 

and note that p > 9 implies p' < O'. 

We prove the upper bound by splitting the sum into two parts. This approach was already 
used in m to prove the result for the Fibonacci cubature rule in W*(T 2 ). Let L n := {d — 
1) log log n and bound the first factor from above by 

0 \ 1 / 0 ' 


( E 2-^Hi|A mjn (-)| ) 

m: \m\i<r„+L n 


V([°d] d ) 


+ 


e\ 1/0' 


va°di d ) 


A. 


( £ 2- 8 '*H.jA mi „(.)| ) 

m: |m|i>r„+L n 

Using p' < 9' and s < 1/6 we use Lemma 13.51 and Lemma 12.31 to bound the first summand by 

< n -V*( r E"(£+ l)‘*- 1 2- fl7( - 1/e) ) W 

e=r n 

< n- s (logn) (d - 1)(1 - s) . 

In the case s = 1/9 this sum is bounded by n~ s (log (loglog n)^ s \ 

To bound the second summand we replace the ^/-norm inside by a £ p >- norm. We obtain for 
s > 1/p again by Lemma 13.51 and Lemma 12.31 the upper bound 


2-S's\m\i 

m: \m\i<r n +L n 


0 ' 

V([0d] d ) 


E 2 -p,s l m li 

m: \m\\>r n +L n 


A. 


m,n 


V([°d] d ) 


w 


< 


n 


Vp( E (£ + l) d_1 2 _p,£ («- 1 /p)^ 

£=r„+L n 
< n- s (logn) (d - 1)(1 - s) . 


W 













































FROLOV’S CUBATURE 


17 


This finally proves Theorem 15.21 


Remark 5.6. It is also possible to study a more general scale of function classes, where the 
mixed smoothness s is replaced by a mixed smoothness vector s = (si,..., Sd). In a way, one 
assumes a bounded mixed smoothness with different derivatives in different directions. Note, 
that this framework is different from the classical “anisotropic” setting in [171 Chapter IV.6] 
and [32], Sect. 4.3], Our definition of the spaces, see Definitions 14.21 &; 14.31 can be adopted 
to this setting by replacing the weights 2 s ' [m ^ e there by 2^ SimiH We denote these 

spaces by F* g . With the standard technique from [131 Lem. C] in combination with our proof 
technique in Subsection 15.11 one can show for 1 < p < oo, 1 < 9 < oo and s > max{l/p, 1 /9} 
the relation 

e(Qn, F B p ,e) ~ n-^Pogn)^- 1 )* 1 - 1 /*) 

where s m i n := min, Sj and v := #{j : Sj = s m in}- Note, that there is no d-dependence in the 
rate anymore. The corresponding results for Besov spaces are stated in [Sj. 


6. Quasi-Banach and limiting cases 

In this section we deal with the remaining cases of the Besov and Triebel-Lizorkin scales, 
that are not treated in the previous sections. We are interested in numerical integration and 
hence only in classes of continuous functions. Besides the quasi-Banach cases min{p, 9} < 1 
with s > 1/p we will also consider the limiting case s = 1/p. In the latter case additional 
assumption are needed to assure continuity. 


6.1. The situation min {p,9} < 1 and s > 1/p. In this section we deal with the classes 
and Fpg with p, 9 < 1, i.e. the quasi-Banach cases. Here, p < 1 affects the asymptotical error 
order negatively, while 9 < 1 does not. The presented lower bounds are given in the upcoming 
Section [71 


Corollary 6.1. Let 1 < p < oo (with p < oo in the F-case), 0 < 9 < 1 and s > 1/p. 

(i) Then 

e(Qn, B 8 p ,e) ~ n~ s . 

(ii) If 's > 1 then 

e(Qn, F; !0 ) X n~ s . 

(in) If 1/p < s < 1 then 

e(Qn, F* p ,o) < n-pogn)^ 1 -). 


Proof. The upper bounds follow from the embeddings A* j for 9 < 1, together with 

Theorems 15.11 & 15.21 


□ 


Corollary 6.2. Let 0 < p < 1, O<0<oo and s > 1/p. 
(i) Then 


e(Q„, B a pt0 ) x n - s + 1 /p- 1 (i og n )( rf - 1 K 1 - 1 / 0 )+ 


(ii) and 


e(Q n , FJ_ # ) x n- +1 /p-i 
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Proof. The stated bounds are a direct consequence of Theorem 15.11 and Theorem 15.21 In fact, 
from Lemma 14.71 we know that we have for p < 1 the embeddings 


T}.s . tW+1 1 /P 

B P ,e ^ **i,0 


and 


■ p,o 


pS+l-l/p As+l-l/p 

*1,1 - B l,l 


□ 


Remark 6.3. In contrast to the case p > 1, there is no dependence on 9 in (ii). In particular, 
there is no effect of “small smoothness”. 

6.2. Limiting cases. In this section we deal with the limiting situation s = 1/p. When 9 < 1 
in the B-case or p < 1 in the F-case we have the continuous embedding A s pd (M. d ) ^A C(R d ). 

Theorem 6.4. (i) Let 0 < p < oo and 0 < 9 < 1. Then 

e(Q n , B$) X n -Vmax{ P ,l}_ 

(ii) Let 0 < p < 1 and 0 < 9 < oo. Then 

e{Qn, F$) X n~\ 

Proof. Due to the embeddings F 1 Jq , B l Jg ^A F]^ = Bj^ if p < 1 it suffices to prove (i) in case 
1 < p < oo, s = 1/p and 9 = 1. It follows the same line as the proof in the case s > 1/p. In 
fact, (15.31) for s = 1/p shows 

|/(/) - QnU) I < n_1/p E 2 |m|l/p ||4 f m * f\\ p < n-Vv \\f\\ 1/p . 

— , p, i 

meN[j: 

|m|i>r„ 

□ 


Unfortunately, our proof techniques do not seem to work for p = 1 in (ii). The following 
result, which is probably not sharp, follows from the embedding in Lemma l4.7I iii). together 
with Theorem I6.4U L We conjecture that e(Q n , F( 0 ) x n _1 (log We leave this 

as an open problem. 

Corollary 6.5. Let 0 < 9 < oo. For any e > 0 there is a constant c e such that 

e(Qn, Fgf?) < c £ n _(1_e) . 

Proof. For an arbitrary p > 1 we have the embedding Fj 0 Bsee Lemma l4.7l iiil. By 
Theorem 16.41 we obtain the rate n -1//p . 

□ 


7. Lower bounds 

The lower bounds that we want to present are valid for arbitrary cubature formulas. For 
this we study the quantity Int n (Frf) from (11.11) for the spaces F^ = A s pd . There are already 
lower bounds for Besov classes in the literature see [45j and the recent works [S3- None 
of those references gives lower bounds for p < 1. We will provide them using an approach 
which is close to the one in [9]. We will use the modern tool of atomic decompositions [58] . 
see Section rrn below, to construct appropriate fooling functions. 
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7.1. Atomic decomposition. We will describe the notion of an atom first. For j £ Nq and 
k £ let Qj k denote the cube with center (2~ 3l mi, ..., 2 ~ 3d m d ) and with sides parallel to the 
coordinate axes of length 2 -J1 ,..., 2~ 3d . For 7 > 0 we denote with 7 Qj )k the cube concentric 
with Qjk with sides also parallel to the axes and length 'y2~ J1 , ,..^2~ 3d . 


Definition 7.1. Let K £ No, L + 1 £ No and 7 > 1. A ii-times differentiable complex-valued 
function aj t k is called ( K , L)-atom centered at Qj.k if 

(i) supp a jik C 7 Q jjk , 

(ii) \D a aj t h{x)\ < 2 a ' 3 , for all a = (or, ... ,a d ) £ Nq, 0 < a* < K , 

(Hi) and there are the coordinate-wise moment conditions 



if i = 1 , 


d, £ = 0, ...,L, ji > 1 . 


The following Proposition is due to Vybfral [58j. Recall, that a p := max{0,1 /p — 1} for 
0 < p < 00 . 


Proposition 7.2. Let 0 < p, 9 < 00 and r £ M. Fix K £ No and L + 1 £ No with 

K > (1 + [rj)+ and L > max{—1, \ a p — rj} . 

If {Xj t k}j,k is a sequence of complex-valued coefficients and {aj t k(x)}j t k a collection of (K,L)~ 
atoms centered at Qj.k then the function 

(7- 1 ) / : = 5Z X hk a 3 .k{x) 

je Ng fee z d 


belongs to g if the right-hand side in 117.211 below is finite. Then it holds 


(7.2) 



< 

r^j 


^ ^2 2l- 7 l 1 ( r_1 / p i 61 



i/e 


Consider a function ip £ C^°(M) with support in [0,1]. Now we define 

a jt k(x 1 ,... ,x d ) := ip(2 jl xi - k±) • ... • <p(2 jd x d - k d ) , x £R d , j £ Nq, k £ Z d 

and observe that aj. k is ( K , —l)-atom for any K £ No centered at Qj.k with 7 = 1 according 
to Definition rm Note, that we do not have moment conditions. According to Proposition [T72] 
we do not need moment conditions if r > cr p . 


7.2. Test functions. Now we are in a position to define test functions of type H7.1D in order 
to prove the Theorem below. By (|7.2p we are able to control the norm || • ||b s • Following (9j, 
Thm. 4.1] we will use test functions of type 

(7.3) a},. := C2-™m-V- 1 V e Jj Jj Ht , 

|j|l=m+l keKj(X n ) 

where D£ := {1, ...,2 ( ' — !},££ N, and Kj(X n ) C Dj := Dj 1 x ... x Dj d depends on the set of 
integration nodes X n := {x 1 , ...,x n } . 

Let A* e fl C'(M d ) denote the function class A* d fl (7(1^) equipped with the norm || • ||a* g ■ 

Theorem 7.3. Let 0 < p, 9 < 00 and s > a p := max{0, 1/p — 1}. 

(i) Then 


Int n (Bp 0 n C(R d )) > n- s+(1/p - 1)+ (logn) (d - 1)(1 - 1/0) + . 
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(U) If 1 < p < oo then 


(in) If 0 < p < 1 then 


Int n (F* 0 n C(R d )) > n ~ s {log „)(«*-i)(i-V«)h 
I nt n (F^nC(M d )) >n- s+1 /v~ l . 


Proof. We only need to prove (i). In fact, (ii) follows from (i) and the embeddings B s 9 e = 
Fg e F* e if 6 > p (Holder’s inequality) and g F s p g if 9 < p (Lemma l4~7|) . The relation 
in (iii) follows from (i) and the embedding B s . s ^ F® where minfp, 9} < 1. 

Let us prove (i): We follow the arguments in m Thm. 4.1]. Let n be given and X n = 
{a; 1 ,..., x n } C [0, l] d be an arbitrary set of n points. Without loss of generality we assume that 
n = 2 m . Since Qj k LI Qj,k' = 0 f° r k k 1 we have for every |j|i = m + 1 a set H)j(X n ) c Bj 
with ffU)j(X n ) > 2^ and X n n B j(X n ) = 0. 

(a) Let p, 9 > 1. We choose the test function (|7.3D with Kj(X n ) := Bj(A' n ). Note, that 
for every k G Kj(X n ) we have that supp aj_ k C [0, l] d . The function g\ e in (17.31) is defined 
via a finite sum of continuous functions with support contained in [0, l] d . Therefore, g\ e is 
continuous and supp g\ 0 C [0, l] d . By Proposition 17.21 and (17.21) we can arrange C > 0 such 
that ||< 7 g 6 >IIb* e < 1 • Clearly, g) e belongs to B® 0 D C(R d ). It is obvious that 


/ dl.o dx 

J [o,i] d 




Of course, a cubature rule admitted in (11.11) that uses the points X n produces a zero output. 
This proves (i) in case p, 9 > 1. 

(b) Let p > 1 and 0 < 1. Let us choose a j G N d with |j|i = m + 1 and define the function 


(7.4) 


9t ■= C 2 - 


/ . a j,k- 

k£j])j (Xn) 


Again, by Proposition 17.21 there is a C > 0 such that \\g 2 ||b s g < 1. Moreover, we have 

f g 2 s dx~ 2~ sm . 


'[ 0 , 1 ]“ 


With the same reasoning as in (a) this proves (i) in case 9 < 1. 

(c) Let p < 1 and 9 > 1. We define the test function 

]T a j>kj , 

|j|l=m+l 

where kj is chosen from Bj(X n ). By (17.21) we find aC>0 such that ||^ 0 ||b s < 1. Computing 

the integral gives 


/ 

J [0,1]' 


g 3 e dx x 2 (- s + 1 /p-i)"* m O-i)a~VO 

1 ’ 


With the same reasoning as in (a) and (b) this proves (i) in case p < 1 and 9 > 1. 

(d) Finally, if 0 < p, 9 < 1 we simply take one single atom gj p := 2~ ms 2 with j G N d 
and k G B j{X n ). Again, there is a C > 0 such that || 5 s iP ||b s 9 < 1 and 


'[o,i] a 


at, P dx 


P ,e 

2 (-s+l/p-l)m 


which finishes the proof. 


□ 
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Remark 7.4. Note that the lower bound in (ii) differs from our upper bounds in the case of 
“small smoothness” p > 6 and 1/p < s < min{l, 1/#}. Hence, to prove optimality of Frolov’s 
cubature formula for each A* () it remains to prove the corresponding lower bound for “small” 
smoothness in Triebel-Lizorkin and Sobolev spaces. This seems to be very delicate and we 
leave it as an open problem. However, we conjecture that our upper bounds are tight. This 
is supported by the fact that for d = 2 the corresponding lower bound was proven for the 
Fibonacci cubature formula, see m Theorem 2.5], which is conjectured to be optimal, cf. 211. 
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